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The di-pion spectrum for the T(nS) — » T(n'S)) transition with n < 4 has the form ^gf ~ (phase space) 



l»7 



- x| , with a; = 



M% n , and AM = M{nS) 



(AM) a -4m; < 5 

the spectrum is shown to reproduce the experimental data for all 3 types of decays: 3 
with r\ « 0.5; 0, and —3, respectively. 
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M(n'S). The parameter 77 is calculated and 
1, 2 -> 1 and 3 -» 2 
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1. The di-pion decays of heavy quarkonia are stud- 
ied for the last three decades, (see pQ, [5] for references 
■and discussion). The first observed process V(2S) — ► 
V(lS)wir, V = ip or T yields a simple di-pion spectrum, 
with the amplitude M ~ (</ 2 — 4m 2 ), and theoretical 
methods have been envisaged (see [2] for a review and 
references), based on PCAC and the multipole gluon 
field expansion (MGFE). However, the 3 — > 1 and 3^2 
di-pion transitions in bottomonium, observed by CLEO 
[3], show two other types of spectra: a double peaked 
spectrum for (3 — » 1) and shallow form for (3 — > 2). Nu- 
merous modifications of MGFE and additional models 
.were suggested [TJH], without, however, a unique phys- 
ical picture for all 3 types of decays. It is a purpose 
of this short note to sketch the general mechanism of 
di-pion transitions, leaving details to the full text in pQ. 

2. Our starting point is the Field Correlator Method 
(FCM) [4], which is based on the use of Gaussian 
(quadratic) field correlators, applicable, in contrast to 
MGFE, also for systems of large size, R A , where 
A < 0.2 fm is the vacuum correlation length 0]. As a re- 
sult all large systems display linear confinement and the 
string tension a defines the dynamics instead of gluon 
condensates as in MGFE. It is important that all de- 
caying X(n) states have R > 0.4 fm > A, and therefore 
should be treated in FCM rather than in MGFE ap- 
proach. 

The pion emission in FCM comes from the light 
quark loop and is decribed by the quark-pion La- 
grangian 5 



with Mb r treated here and in pQ as a fitting parame- 
ter, found from the pionless decay, and U(x) is a chiral 
matrix 



U (x) = exp ij 5 
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Note, that U (x) describes creation of any number of 
ir(K, rj) through small dimensional factor /„. = 93 MeV 
in the denominator. This ensures strong interaction and 
possible multipion-quarkonuim resonances. 

The QQ Green's function with light quark loop in- 
side and with two possible ways of emission is shown 
in Fig. 1(a) and (b). We shall keep here notations "a" 
and "b" for the amplitudes with subsequent one-pion 
and zero-pion-two-pion emissions, respectively. 

Using notations n and n' for X(nS) and X(n'S) 
states, and 712,713 the Qq and Qq states, respectively, 
one can write for the amplitudes a and b, 
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Sqm = —i d xtp(x)M br U (x)ip(x) 
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Here 7 = J,J^,J^ are overlap integrals of 

heavy quarkonia wave functions i^qq and the product 



1 



2 



of Qq and Qq wave functions tpQq' "^Qq > multiplied by 
the exponentials of free relative motion exp(zpr) with 
momentum p, and the pion plane waves with momenta 
ki and k2. 

The resulting decay probability is 



dw ^ 
dq 



n') = /i 2 JW^x)d cos 9\a-b\ 2 

(5) 
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where y? = (AM) 2 — 4m 2 , x = 
(fci + &2) 2 and 9 is the angle of emitted ir + with initial 
X(nS) direction. It is important, that a and 6 depend 
differently on q, 9: in b, the vectors ki,k2 enter as a 
sum K = ki + k 2 , while in a one has a product of de- 
creasing functions of |ki| and | lt2 ] - E.g. m the Simple 
Harmonic Oscillator (SHO) basis one can write a form, 
exactly satisfying the Adler zero condition: 
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T nn ' (^1 + U 2 ) + UJ\ + UJ 2 



Here 



and /?2 is the SHO parameter 



for the B, B* wave functions. Fitted to known r.m.s. 
radius one has /3 2 = 0.5 GeV. Now one can see, that 
k 2 + k 2 . = a(q) + ^y(q) cos 2 9 is a weak but increasing 
function of q, while K 2 is strongly decreasing function 
of q. Expanding a, b in powers of x at the threshold, 
one arrives at the expression: 
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dcos9y / x(l — x)b 2 h \rj— x\ 



(7) 

Here n is a real parameter below BB threshold. We 
have calculated the values of n for Y(nS) — > T '(n' S)tttt , 
using the SHO basis with parameters fitted to all r.m.s. 
radii of states and finding 77 from ([6]), where r nn i were 

found from The intermediate states in all three 

transitions (3 — > 1), (2 — > 1), and (3 — > 2) were chosen 
as (BB* +B*B) for the amplitude a, and (BB) for the 
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amplitude b. We have checked that the Adler zero con- 
dition is fulfilled when the full set of intermediate states 
is involved, and imposed the Adler zero condition in our 
case of restricted set of intermediate states, i.e. we have 
used the form ([6]) with a(k\ = u>i = 0, k 2 ) = b(k\ = 
u>i = 0, fc 2 ). 

As a result we have obtained the following values 
of r\: n(3 1) 0.39, rj(2 -> 1) = 0.05 and 
77(3 — ► 2) = —3.2. On the other hand we have inde- 
pendently fitted the three di-pion spectra measured by 
CLEO using the form 

dw(n — > n') 



dq 



const. V 'x(l - x)\r)fit - x\ 2 (8) 



with constant and rjtit as free parameters and found 
that nfit(3 — > 1) w 0.5, 

r]fit(2 — > 1) = 0, ??/it(3 —> 2) ps —3. Results of this 
fitting are shown in Fig. 2 

3. The fitting curves in Fig. 2 appear to reproduce 
the main features of di-pion spectra in all three cases. 
The deep well in the (3 — > 1) double-peaked fitting curve 
is partially filled since n actually depends on cos9 and 
in integration J dcos9\n(9) — x\ 2 one has nonzero result 
for all x. 

It is interesting that the total yield of nn, calcu- 
lated for ip(2S) -> J/ipmr with M br = 1 GeV (fitted 
to ^(3770) — > DD), yields the value of which is 3 
times smaller than experiment. This is reasonable, since 
Ttttt is proportional to the fourth power of the overlap 
integral and the SHO wave function is certainly a very 
rough approximation, nevertheless, the resulting is 
in the correct ballpark. 

The di-pion spectrum for ip(2S) — > J/ipirir also 
corresponds to 77 w 0, while those of T(4S f ) — > 
T(n'S) nn with n' = 1, 2 have J7/ it (4S -> 25) = 0.30 and 
nfit(AS — » 25) = 0.61 respectively. The calculation of 
?7(4S' — > n'S"), however, should be done not with SHO, 
but with realistic wave functions, which is in progress. 

The processes T(55) -> T(n'S)-mT with ri = 1,2 
were observed recently [6] with di-pion spectra of possi- 
bly three-peaked form: also the total di-pion yield Y^* is 
thousand times larger, than for n — 4S, 3S. Both facts 
can be derived from our expressions, Eqs. (|4l6p . if one 
takes into account that the mass of T(55) = 10.86 GeV 
is above all three thresholds BB, BB* , B* B* and the 
appearing imaginary part is very large. This analysis is 
now in progress. 

Finally, one can compare the total yield of K + K~ 
to that of 7r + 7r - in T(55 — > IS) transition. Indeed, the 
only difference is in /if = (AM) 2 — 4m 2 , i = tt,K. 
From ([7]) one can derive that Ti — f ^dq ~ nj, 



and 



r x ,(5S^is) 
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0.104, which roughly 
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agrees with experimental data (an additional 50% re- 
duction of this ratio follows if one use fx — 0.112 
GcV^ U = 0.093 GeV). 

Summarizing, we have suggested a new nonpertur- 
bative approach which describes all transitions X(n) — > 
X{n')im(KK) for n = 2, 3 and n' = 1, 2 without fitting 
parameters. It is argued that the standard MGFE ap- 
proach to di-pion transitions in heavy quarkonia cannot 
be applied to large size quarkonia. 
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Fig. 2 The experimental data for T(nS) — ► T(n'S)irir 
from [3] versus theoretical spectrum, Eq([7]) with r\ = 
0.55; 0; —2.7 (top to bottom, dashed line). 



